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I .  INTRODUCTION 

There  are  presently  only  two  basic  analytical  methods  applicable  to 
scattering  from  random  rough  interfaces  --  the  Rayleigh  or  small  wave -height 
approximation  and  the  Kirchhoff  or  physical  optics  approximation.  Both  of  these 
methods  are  valid  only  for  a  limited  range  of  parameters.  Composite -roughness 
theory,  which  combines  these  two  approaches  in  their  respective  regions  of 
validity,  also  has  a  defect;  its  predictions  are  dependent  on  an  arbitrary 
parameter . 

The  Rayleigh  approximation  arises  from  assuming  that  the  scattered  field  can 
be  written  in  terms  of  waves  propagating  away  from  the  surface,  even  in  the 
hollows  between  adjacent  peaks  on  the  surface.  The  field  on  the  bounding 
interface  is  not  explicitly  computed.  The  Kirchhoff  approximation,  on  the  other 
hand,  arises  from  retaining  only  the  lowest  order  term  in  the  solution  of  an 
integral  equation  of  the  second  kind  for  the  surface  field.  This  approximate 
surface  field  is  then  used  to  compute  the  scattered  field  via  the  Helmholtz 
integral . 

In  addition  to  the  integral  equation  of  Mane1  and  Meecham2,  which  is 
traditionally  used  when  discussing  the  Kirchhoff  approximation  for  scattering  from 
pressure  release  surfaces,  there  exist  other  surface  field  formulations:  namely, 
Uretsky's3  integral  equation  of  the  first  kind  and  a  variant  of  the  Maue-Meecham 
equation  introduced  by  Dashen4.  In  this  article,  the  solution  of  these  three 
different  integral  equations  is  obtained  in  the  form  of  a  series  by  first 
expanding  the  kernel  of  the  integral  equation  in  a  small  parameter  and  then 
applying  the  method  of  successive  approximations. 
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Section  II  introduces  the  problem  addressed,  scattering  from  a  one¬ 
dimensional  periodic  pressure  release  surface  and  reviews  the  result  obtained  when 
the  Rayleigh  hypothesis  is  invoked.  In  Section  III,  the  three  integral  equations 
are  solved  and  expressions  are  obtained  for  the  reflection  coefficients. 

Uretsky's  equation  yields  a  result  identical  to  the  Rayleigh  series  to  third 
order.  The  two  integral  equations  of  the  second  kind,  however,  yield  results  that 
differ  from  each  other  and  the  Rayleigh  series  at  third  order.  Only  one  of  these 
series,  that  based  on  Uretsky's  integral  equation,  has  the  necessary  attributes  to 
represent  a  meaningful  physical  solution.  The  reflection  coefficients  arising 
from  the  two  integral  equations  of  the  second  kind  are  not  reciprocal  and  also  are 
unbounded  when  the  grazing  angle  associated  with  the  reflected  order  vanishes. 
Uretsky's  equation  may  also  be  cast  into  an  integral  equation  of  the  second  kind 
and  the  convergence  of  its  solution  demonstrated  by  a  method  due  to  Urusovskii5. 
However,  if  the  kernel  of  this  integral  equation  is  expanded,  this  proof  is 
inapplicable.  This  is  also  true  for  the  other  two  integral  equations  leading  to 
the  disparate  results. 

Section  IV  summarizes  these  findings  with  the  conclusion  that  Uretsky's 
integral  equation  is  the  preferred  starting  point  for  the  development  of  improved 
low-frequency  scattering  formulations. 

II.  THE  SMALL  WAVE  HEIGHT  APPROXIMATION 

Figure  1  depicts  the  problem  of  interest  --  scattering  of  an  incident  plane 
wave  of  wavelength  A  -  2rr/k 

Pinc(x,z)  “  exPtik(Q0x  -  T0z)].  (1) 

by  a  one-dimensional  periodic  pressure  release  surface.  Our  nomenclature  closely 
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follows  Holford's6  with  an  incident  grazing  angle  ^  so  that  -  cos^^  and 

7q  -  sin^Q.  Far  from  the  surface,  the  total  field  may  be  written 

p(x,z)  -  p.  (x,z)  +  S  R  exp[ik(a  x  +  7  z ) 1 ,  (2) 

^  rtnc  nn  n  n 

where  R  is  the  reflection  coefficient  for  the  nC^  order  and 
n 

<*n  -  aQ  +  (nK/k)  , 


a  -  ^  : 

‘  <4-  1>1/2  : 


where  the  fundamental  surface  wavelength  A  -  2x/K.  Note  that  in  Eq.  (2)  when  the 
limits  on  the  summation  are  omitted,  the  limits  are  understood  to  be  -<*>,  +<*>.  This 
same  simplification  will  also  be  used  for  infinite  integrals. 

Rayleigh  assumed  that  Eq .  (2)  held  not  only  in  the  far  field  but  also  on  the 
boundary  z  -  f(x).  Equating  p[x,  f(x)]  -  0,  then  yields 

exp[ -ik7QC(x) ]  -  S  R^fiKnx  +  iky^x)].  (4) 

To  obtain  a  solution  to  this  equation  R^  is  expanded  in  a  small  parameter 
e  -  0  (kh) 

R  -  R(0)+  e  R(1)+  € 2  R(2)+ .  (5) 

n  n  n  n 

The  exponentials  in  Eq.  (4)  with  arguments  -  0  (c)  are  also  expanded,  permitting 
an  iterative  solution  for  the 

n 

R(0)-  -6  .. 


2ikW 


2k  7n  S  7  h  h 
Om  m  n-m  m 


(ik37  /2)  S  h  h  h  l-(7^+  ih/'i  +  12+  l2-  7  ]. 

U  m,p  n-p  p-m  m  n  u  m  p  n>  p 
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with  h 


0JA  5"(x)  exp(-iKmx). 


The  identity 


n2Sh  h  h  =3  S  h  h  h  (p2-  m2),  (7) 

n-p  p-m  m  m,p  n-p  p-m  m  K 

m,  p  r  r  r 

which  can  be  verified  by  interchanging  the  indices  has  been  used  to  cast  the  third 
order  reflection  coefficient  into  a  reciprocal  form. 


III.  SURFACE  INTEGRAL  FORMULATIONS 

A.  The  Maue -Meecham  Integral  Equation 

Maue1  and  Meecham2  derived  an  integral  equation  of  the  second  kind  for 
the  normal  derivative  of  the  pressure  on  the  boundary 


0(x)  -  20(x)  -  /  0(x')  K  (x',x)dx', 


where  0<x)  -  ~  ( -  fz  +  f' (x)  ^  ), 


V>(x)  =  -[70+  a0f'(x)]  Pincfx>C(x)  ]  , 
ik  H{1)(kp) 

K^x'.x)  -  J - - - Ef(x')  -  f(x)  -  (x'-  x)f'(x)], 

and  p  -  {(x'-  x)2  +  (f(x')  -  r(x)]2)1^2  . 

The  Kirchhoff  approximation  results  if  only  the  first  term  on  the  right  hand  side 
of  Eq .  (8)  is  retained. 

On  applying  the  operator  ~  exp(-ikomx;  to  both  sides  of  Eq.  (8),  one 


obtains  the  set  of  equations 


0  +  E  V(1)  0  =  20 

m  m ,  n  n  m 

n 


where 
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(1)  ik  rA  dx 

m ,  n  2  0^  A 


v\-/  _  iii  ii*  J  dr  (r,x)  expfika^r  +  iK(n-m)x], 


(10) 


*m  "  0-fA  f  ^(X)  exP('ik%x). 

and  r  -  x ' -  x . 

To  obtain  tue  desired  series  solution,  the  Hankel  function  appearing  in  the  kernel 
of  the  integral  equation  is  expanded  in  the  series7 

(Ur)'u/2  H^1)  [  z  ( 1+r ) 1/2  ]  =  S  ^  (“)  (z)  ,  (11) 

m=0 

which  is  absolutely  convergent  for 
I  r  |  -  [f(r+x)  -  f (x) ] 2/ r 2  <  1. 

The  integrations  indicated  in  Eq .  (10)  may  then  be  performed  (see  the  Appendix)  to 
obtain 

V(1)  -  i  h  [k  (7  -7  )  +  («n-n)  K  o  /7  ] 
m,n  m-n  m  n  n  'n 


+  ( ik  /2)  2  h  h  h  {  K[(m-q)  c*  7 

m-q  q-p  p-n  q  q 

q  1  P 

+  (p-n)an7n  -  2(q-p)ap7p]  -(k/3)  (7^-  7^+  37^-  37^) > 

+ .  (12) 

It  is  of  interest  to  note  that  the  leading  terra  in  this  expansion  can  be  written 
in  the  form 


JA  7*  exp[  iK(n-m)xK"  (x)  , 


q  (a  +  a  )  , 

1  .  n  m  n  1 

k  ,  .2  +  ( 7  +  7  )  °J  A 

7  (7  +  7_)  m  n 

n  m  n 


(13) 


from  which  it  is  evident  that  this  is  a  curvature  correction. 
If  is  expanded  in  powers  of  « 
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and  the  exponentials  of  argument  -  0(e)  are  also  expanded,  Eq.  (9)  may  be  solved 
for  the  \£>^  ^ 


2  70  Sn,0 


2itc  7(]ynhn 


k27_  Eh  h  (2  7  7  -  72) 

0  n-m  m  m  n  n 

m 


0^-  ik37n  2  n  h  h  [7n(72  -  72/3  -  27  7  ) 
n  U  n-m  m-p  p  n  m  n  m  p 

m,p 


+  7  (l  +  2a  +  a  a  -2a  a  -2a  a)], 
p  p  p  n  nm  mp 


The  identities 


Eh  h  h  (p)  -  *  Eh  h  h  (m) 

n-m  m-p  p  2  n-m  m-p  p 

m,p  m,p 

=  J  Eh  h  h  (n) 
j  n-m  m-p  p 

m,p 

are  used  in  obtaining  Eq .  (14). 


With  known,  the  reflection  coefficients  are  found  using 
Rm  ”  2t  Z  ^n  f  eXp[iK(n’m)x  'ikV(x)]- 


m  n 


The  R  are  expanded  as  in  Eq .  (15)  and  the  exponential  factor  with  argument 
m 

-  0  (e)  is  likewise  expanded.  For  0  <  j  <2,  the  ^  are  identical  to  the 
Rayleigh  reflection  coefficients.  For  j  -  3,  the  result  differs 

ik^7 

f  3  7  (3)  'o 

R*1  '  (Rayleigh)  +  - -  E.  h  h  .  h. 

n  n  7  6  7n  P.J  n-p  p-j  j 


X(VV  anQp‘2)  ’7j  (Qn  +  QnQj  ’2)  +  (7n  *  70  +  a070Qn)/3]  ’ 


(16) 
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The  additional  terra  in  Eq .  (16)  is  not  reciprocal  and  also  becomes  unbounded 
vanishes . 

B .  Dashen's  Formulation 

Dashen4  introduced  the  integral  equation 


V>(x)  -  -2yn  (JJ)2  P.  [x,f(x)]  -  (tV  J  0(x')K,  (x'  , x)dy '  , 


ds  „  2 


'0  dx  inc 


dx 


(17) 


ik  HJ1)<k'> 

where  K2(x' ,x)  -  — 


and 


O2  - 1  +  ^'<x> 


p 

,2 


:r(x')  -c (x) ] , 


Performing  the  same  operations  as  in  subsection  A,  the  set  of  equations 


4>  +  E  v(2)  -  H  , 

m  m,n  n  m 

n 


(18) 


is  obtained.  In  this  case 

0, 


m  0 


JA  exp(-iKmx)  (^) 2  exp  [  -  ik7nf  (x)  ]  , 


'0s 


and 


Vi2i  ‘  f  oT  f~  <K>"  I  dr  K2(,’x>  e,tpllhV  +  iKCn-n,)^] 


ik  rA  dx  ,dss2 


(2) 

With  the  use  of  Eq.  (11),  ^  may  be  expanded 


V 


(2)  -  ik(7  -  7  )  h  -  ikK2  S  h  .  h.  h  (j-p) 
m,n  m  n  m-n  .  m-j  j-p  p-n 

J,P  J  J 


....  .  ik  ,  ,  ,  .3  3-3  -,3. 

x(m-j)  (7  -  7  )  +  ~  S  h  .  h.  h  (7  -  7  +  37.  -37  ) 

J  p  m  6  .  m-j  j-p  p-n  'n  'm  j  p 


The  first  term  in  Eq .  (19)  may  be  rewritten  as 


(19) 


f 77 n)  0/A  f4  r (x)  exp[iKx(n'm)1, 

v  rm  rn) 


as  7 

n 
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and  is  a  slope  correction. 

The  surface  field  transform  0  can  again  be  expanded  in  powers  of  e.  For 
Dashen's  equation,  the  result  differs  at  third  order  from  Eq.  (14).  In  this  case, 


0. 


(3) 


n 


y  y  3 

-2ik\  E  h  h  h  {-7-  +  Hj- 
0  n-m  m-p  p  6  3 

m ,  p 

2  2  2  2  2 
+  7o  (7p  ‘  7m  +  7n/3  '  V3)  +  7n7m7p'  7n7m/2 

+  (7p/3)  (apaQ  -  aan  +  a^-l)] 


(20) 


The  corresponding  reflection  coefficient 

ik3-y 

„(3)  0(3),„  ,  .  _  .  70 


R'  '(Rayleigh)  +  ~  E  h  h  h 

n  n  y  0  2y  n-m  m-p  p 

n  m,p 


x[2y  7n(7r,'7  )/3  -  4y  7^/3  -  7^7  -  73 

n  0  0  n  p  n  n  m  p 

-2-y  (a  a  -  a  a  +  a  a_-l)/3],  (21) 

p  p  0  p  n  n  0  '  J 

differs  from  that  of  subsection  A,  but,  nevertheless,  suffers  from  the  same 
deficiencies . 

C .  Uretskv's  Integral  Equation 

Uretsky's3  integral  equation  is  obtained  by  evaluating  the  Helmholtz 
integral  on  the  boundary 

J  0(x')  K0(x,x)dx'  =  -2pinc[x,r(x)] .  (22) 


where  KQ(x' ,x)  -  (k/2)  Hg^(kp). 
equations 


E 

n 


V<0)  0  - 

m ,  n  n 


This  integral  equation  reduces  to  the  system  of 

(23) 


where 


0m  =  (/  f  e*P('ikVO  Pinc[x,r(x)]( 
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m ,  n  0 


,/A  A“  f  drKQ(r,x)  expfik a^r  +  iK(n-m)x] 


Again,  using  Eq .  (11),  may  be  written  as  a  series 

V(0)  =  6  /7  -  k2  E  h  h  <7  +7  -2y  )  +  ...  .  (26) 

m.n  m,n  m  m-p  p-n  m  n  p 

The  second  term  in  Eq .  (26)  can  be  manipulated  to  show  that  it  depends  on  both  the 
square  of  the  surface  slope  and  the  product  of  the  curvature  and  height. 

The  surface  field  transform  differs  at  third  order  from  that  of  subsection  A 

0(3)»  ik37n7  Z  h  h  h  [-7 i/3  +  77  +  7  2  -27  7  ] ,  (27) 

n  0  n  n-p  p-m  m  1  Cr  'm'n  m  'p'mJ’ 

m ,  p  r  r 

and  the  reflection  coefficients  are  identical  to  the  Rayleigh  reflection 

coefficients.  In  a  related  study8,  it  is  shown,  for  the  case  of  a  pressure 

release  sinusoidal  surface,  that  che  reflection  coefficients  resulting  from  the 

expansion  of  Uretsky's  equation  are  identical  to  the  Rayleigh  result  at  every 

order  in  the  expansion. 

With  Eq .  (24),  Eq .  (23)  can  be  rewritten 

</>  +  Z  V  V>  -  - 2y  0  ,  (28) 

m  m ,  n  n  mm 

n 

where 


V  -  7  V  -6  .  (29) 

m,n  m  m.n  m,n 

For  the  solution  of  Eq .  (28)  to  converge  as  the  number  of  equations  retained 


becomes  large,  ^  must  be  a  completely  continuous  operator  in  the  space  of 
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square  summable  sequences,  and  the  right  hand  side  of  this  equation  must  be  an 
element  of  this  space,  or9 


2 

m ,  n 


I  27 


0 

m  m 


2 


<  00 , 


(30) 


and 


2  IV  I  <  °°- 
m,n 

m,  n 


(31) 


From  Eq .  (29),  this  latter  requirement  will  be  satisfied  if 


2 

m ,  n 


r  v^0) | 2<  ®. 
m  m ,  n 


(32) 


Urusovskii5  has  shown  that  7  V  ^  is  bounded.  Thus  with  m-n+p ,  Eq .  (32)  becomes 

m  m, n  r  n 


2  1 7  V<°>  |2 

n+p  n+p , n  < 

n ,  p 

f n f , |p|>  N  <  ® 


(33) 


With  the  assumption  that  $"(x)  is  four  times  differentiable,  Eq.  (23)  may  be 
integrated  by  parts,  twice  with  respect  to  both  x  and  r,  whence 


1 7  V(0)  |<  ~  as  |  p  |  ,  |  n  I-  oo, 

n+P  n+P’n  | p | 2 | n | 2 


(34) 


where  M  <  This  is  sufficient  to  satisfy  Eq.  (33).  Two  integrations  by  parts 
of  Eq .  (34)  are  sufficient  to  demonstrate  that  Eq .  (30)  bolds.  Thus, 
Kantorovitch's  method  of  reduction9  is  applicable  to  Uretsky's  formulation  as  well 
as  the  two  integral  equations  of  the  second  kind.  Note,  however,  that  this  proof 
is  inapplicable  to  the  expansions  of  ^  given  in  Eqs .  (12),  (19),  and  (26). 

D .  Monostatic  Backscatter ine  Strength 

The  expressions  given  for  the  reflection  coefficients  are  of  sufficient 
complexity  that  it  is  difficult  to  ascertain  their  dependence  on  various 
parameters.  If  one  takes  ensemble  averages,  the  resulting  expressions  are 
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simpler.  Let  us  first  consider  the  quantity 

Q  .  -  <  R(1VaL  R<2>  R*<2>>  +  <  R^R*1**,  .  (35) 

n,J  n  j  nj  n  j  n  j 

which  is  needed  to  compute  the  backscattering  strength. 

For  a  Gaussian  distribution  of  surface  heights 

<h  h  h  h  >  -  W(Kp)  W(Kr)  S  S 
p  q  r  s  r  P>’Q  r , -s 

+  W(Kp)  W(Kq)  65+  W(Kp)  W(Kq)  6  6 

P .  -  r  q ,  -  s  r  p,-s  q  ,  -  r 

where  W(Kp)  is  the  rough  surface  wave  number  spectrum  normalized  such  that 

<  C2(x)  >  -  S  W(Kp) . 

P 

The  first  term  on  the  right  hand  side  of  Eq.  (35)  is  the  same  for  all  of  the 
formulations 


QnJ  ’  4k2’0S„.jIU<2k“0>  *  lk<“„  -V 

J  J  m 


x  W  (2ka  +  Km)  [77*  +  Re<T  7*  )]),  (36) 

u  m  m  m  -m 

for  monostatic  backscatter,  where  Re(x)  denotes  the  real  part  of  x.  For  the 
Rayleigh  and  Uretsky  formulations,  the  second  term  on  the  right  hand  side  of 
Eq.  (35)  is 

qnj  -  8k^0Sn.j  “<2kV  Op^oWV1' 

For  the  Maue-Meecham  and  Dashen  formulations,  respectively,  one  obtains 

«nj  -  Qnj  -  16o2ik4Vn  J  W(2to0) 

x  Re [ S  W(Kp)  (70-7  )] , 

P  P 

and 
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Qnj  ‘  Qnj  <R“yl^el')  -  ^Vn.j  «'2kV 


x  Re { 2  W(Kp)  [27p(70-7p)  +  8a0(7pQp-7pa0-70a0)/3] 


The  scattering  strength  is  then  given  by 


S 


k 

K 


<}* 


q<2>>. 

n>J 


Figure  2  compares  the  backscattering  strength  predicted  by  the  different 
formulations.  The  wave  number  spectrum  assumed  for  this  computation  is 


~  W(Kp) 


,01 

(Kp)3 


.05/m  <  Kp  <  100/m. 


IV.  CONCLUSIONS 
i 

The  key  to  solving  a  general  scattering  problem  is  first  to  obtain  a  solution 
for  the  relevant  field  parameter  evaluated  on  the  scattering  surface.  For  all  but 
very  simple  cases,  which  can  be  treated  numerically,  this  entails  an 
approximation  which  thence  limits  the  region  of  applicability  of  the  result 
obtained . 

This  article  has  considered  three  different  formulations  of  the  boundary 
integral  equation  for  plane  wave  incidence  on  a  periodic  pressure  release  surface. 
In  the  low-frequency  limit,  the  three  equations  yield  different  results.  Only 
Uretsky's  integral  equation  yields  physically  acceptable  reflection  coefficients 
and  these  are  identical  to  those  obtained  using  the  Rayleigh  or  small  wave  height 
approximation . 

Whereas  the  Rayleigh  approximation  is  in  good  agreement  with  exact  solutions 

C _ — 

for  scattering  at  low  frequencies  and  grazing  angles10,  when  the  frequency  or 
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grazing  angle  is  significantly  increased,  this  approximation  yields  scattering 
strength  predictions  considerably  higher  than  those  experimentally  observed.  It 
appears  that  Uretsky's  integral  equation  is  the  only  viable  starting  point  for 
developing  a  general  scattering  theory  applicable  over  a  greater  range  of 
parameters.  '  ^  _____ 


ACKNOWLEDGEMENT 

This  research  was  performed  at  the  Applied  Research  Laboratory  of  The 
Pennsylvania  State  University  with  the  support  of  the  Office  of  Naval  Research. 


15 


2  March  1990 
STM : lzh 


Appendix 


All  of  the  integrals  considered  are  evaluated  using  the  following 
prescription.  The  Hankel  functions  are  first  expressed  in  the  form 
i  j  — & 


H(1)(Z) 

V 


iz  <k2-  t>2)1/2  .  .  -l.(k2-  o2)1/2. 

2  2  1/2  expl k  -  l“  sln  1  k 


(k  -  p  ) 


(A- 1) 


Terms  of  the  type  [C(r+x)  -  f(x)]  are  expressed  as 

da  E  h^  (iKmr)  exp [ iKm(x+ar ) ] ,  (A-2) 

m 

and  those  of  the  form  [^(r+x)  -  f (x)  -r{"(x)]  as 

-f1  "d/9  S  h  (iKmr)2  exp(  iKm[x  +  (l-^)r])  .  (A-3) 

U  m 

m 

By  assigning  k  a  small  positive  imaginary  part  5,  the  order  of  integrations  can  be 
reversed  and  the  integration  over  r  performed.  The  integration  over  p  in 
Eq.  ( A- 1 )  may  then  be  performed  by  closing  the  contour  in  the  upper  half  plane,  at 
which  point  one  can  set  6=0.  The  remaining  integrations  are  straight  forward. 
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Figure  2.  Monostatic  backscatter ing  strength  predicted  by 
expansion  of  the  three  integral  formalisms. 
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